Statistics of orbital entanglement production in quantum-chaotic dots 
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The production of orbitally entangled electrons in quantum-chaotic dots is investigated from a 
statistical point of view. The degree of entanglement is quantified through the concurrence and the 
entanglement of formation. We calculate the complete statistical distributions of the entanglement 
measures by using random matrix theory. Simple analytical expressions are provided for the con- 
currence distributions. We identify clear signatures of time-reversal invariance in the production of 
entanglement at the level of the entanglement-measure distributions, such as the ability of producing 
maximally entangled (Bell) states, which passed unnoticed in previous works where only the first 
two moments of the distributions were studied. 
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In the early times of quantum mechanics, the predic- 
tion of entangled quantum states subject to nonclassi- 
cal correlations was identified as a possible sign of in- 
completeness of the theoryji Despite such original objec- 
tions, and after many years of debate during which deci- 
sive experimental evidence was accumulated, the reality 
of quantum entanglement is nowadays widely accepted. 
Besides its fundamental interest, modern research on en- 
tanglement has been strongly motivated by its possi- 
ble application as a resource for quantum information 
processing^ As a consequence, an intensive research ac- 
tivity on the production, manipulation, and detection of 
quantum entanglement in a variety of physical systems 
is reflected in an extended literature. 

Several efforts have been devoted to the study of the 
entanglement in electronic systems^ with a view on 
solid-state applications that could lead to a major tech- 
nological breakthrough in the field of nanoelectronics. 
As far as the production of entangled electrons is con- 
cerned, several proposals based on different interacting 5 
and noninteracting^i 7 -^ electron mechanisms already ex- 
ist. Recently, Beenakker et. al^ proposed a ballistic 
quantum dot as an orbital entangler for pairs of nonin- 
teracting electrons^; by assuming that the entangler is 
a chaotic quantum dot, they calculated the average and 
variance of the concurrence (a standard measure of two- 
qubit entanglement) and found that these two moments 
are practically unaffected by the breaking of time- reversal 
invariance (TRI). This fact is in contrast to other trans- 
port properties of ballistic quantum dots such as the con- 
ductance, where TRI yields significant weak-localization 
corrections. More recently, it was shown£ that signatures 
of TRI corrections can arise in the average concurrence 
once spin-orbit interaction is considered. The concur- 
rence fluctuations, however, can be very large (of the or- 
der of the average concurrence).— This indicates that the 
first two moments of the concurrence are insufficient for 
an accurate statistical description of the entanglement 
production in a quantum dot. 



Here we calculate the complete distribution of the de- 
gree of orbital entanglement in terms of the concurrence 
C and the entanglement of formation £. We derive sim- 
ple analytical expressions for the distribution of C for the 
cases when TRI is preserved and broken, while the distri- 
butions of £ are calculated via a relation between C and 
£ . We find that the degree of produced entanglement is 
distributed quite differently depending on whether TRI is 
present or not. In particular, we show that maximally en- 
tangled states are produced only when TRI is preserved. 
Clear signatures of TRI correlations are, however, absent 
in the first two moments of the distributions, which is in 
agreement with previous results* 7 - We verify our theoreti- 
cal results for the concurrence distributions by numerical 
simulations of chaotic quantum dots. 

The setup of the orbital entangler as proposed in Ref. 
is sketched in Fig. Q] It consists of a quantum dot with 
two attached single-channel leads at the left and right. 
Each lead is connected to an electron reservoir. The ap- 
plication of a small bias voltage between reservoirs give 
rise to a coherent current traversing the dot from left to 
right. Exchange correlations due to scattering within the 
dot lead to entanglement between transmitted (to the 
right) and reflected (to the left) electrons, as we show in 
the following. 

We start by considering a separable two-particle state 
incoming from the left reservoir in Fig. [T] 

|* in )=aI4|0), (1) 

where a\ creates an incoming electron in the lead i = 1, 
and 2 and |0) denotes the Fermi sea at zero temperature. 
We disregard spin degeneracy for simplicity (equivalently, 
one can consider incoming electrons with their spins po- 
larized along the same direction). The outgoing state is 
a coherent superposition of orbital channels determined 
by the single-particle scattering matrix 5, such that 

\^out)=Y, S ^Sn2blbi\0). (2) 
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FIG. 1: The quantum dot orbital entangler. The left and 
right leads are attached to chemical potentials pl and (ir, 
respectively. As explained in the text, an electron leaving the 
dot to the left (right) side can escape through the leads 1 or 
2 (3 or 4), defining a two-level quantum system. 



Here, &j creates an outgoing electron in lead j, while the 
matrix element Sji describes the scattering of an incident 
state from the lead i(= 1 and 2) to the lead j(= 1, . . . , 4). 
Notice that the terms with m = n in Eq. ^ vanish for 
the sake of fermionic statistics. The state l^out) can be 
split into three components representing sectors of the 
Fock space with different local particle number at the 
left (til) and right (tir) of the dot such that the total 
particle number remains constant (til + «r = 2): 



Y, K,n R ) = |2,0> + |0,2> + |1,1) 



with 



|2,0) 
|0,2> 
11,1) 



(S 31 S i2 - S 41 S 32 )blb{\0), 
(S 11 S 22 -S 2 iS 12 )b\bl\0), 

^^(SplS q2 — S q l S p2 )blb\\0), 



(3) 



(4) 
(5) 
(G) 



where p — 1, 2 and q = 3, 4 in Eq. ©. We intend to char- 
acterize the degree of orbital entanglement between the 
left and right outgoing channels. In practice, this could 
be addressed by studying current correlations in Bell- 
like measurements involving local transformations at the 
left and right sides of the dot^ L 10 ' 11 Such transforma- 
tions acting on the state |^ ut) conserve the local particle 
number, forbidding a mixing between the different terms 
of Eq. ([3]). Under these conditions, the total accessible 
entanglement^ is determined by a sum of distinct contri- 
butions from each individual sector of the Fock space. In 
our case, the sector contributing to the orbital entangle- 
ment is the one with equal occupancy at both sides of the 
dot (i.e., |1, 1)) since |2, 0) and |0, 2) are clearly separable 
in terms of the bipartition left-right £L£ An electron leav- 
ing the quantum dot to the left side can choose between 
the leads 1 and 2 for escaping (see Fig. [I}. This defines 
a two-level quantum system or qubit. The same happens 
with an electron escaping to the right side through leads 
3 and 4. This means that the component |1,1) in Eq. 
(|rj|) describes (up to a normalization factor) a two-qubit 
entangled state. 



A widely used measure for quantifying two-qubit en- 
tanglement is the concurrence C. It is defined as~2- 



C(p) = max{0, Ai — A2 — A3 — A4}, 



(7) 



where the Aj's are the eigenvalues (in decreasing order) 
of the matrix pp, with p being a 4 x 4 two-qubit den- 
sity matrix (p = |1, 1) (1, 1|/(1, 1|1, 1) in our case) and 
p = (tj y ® <? y )p*(o'y <B> cr y ), where <j v the second Pauli ma- 
trix. The concurrence varies from to 1. The case C = 
corresponds to separable nonentangled states, while max- 
imally entangled (Bell) states own C = 1. Those states 
with a < C < 1 are non-separable partly entangled 
states. For our quantum dot entangler in Fig. 1, a fi- 
nite value of C would guarantee that the left and right 
outgoing channels are orbitally entangled. 

The concurrence has the advantage of being directly re- 
lated to the entanglement of formation £, which is one of 
the most accepted measures of entanglement. Physically, 
£ quantifies the cost, in terms of Bell states, to prepare 
a given (pure) state. The entanglement of formation and 
the concurrence are related through^ii^ 



£(C) 



where 



h(x) = — x\og 2 {x) — (1 — x) log 2 (l — x). 



(8) 



(9) 



On the one hand, C and £ can be written in terms of the 
scattering matrix S through the transmission eigenvalues 
Ti and t 2 as&£ 



C = 



2^r 1 (l-r 1 )r 2 (l-r 2 ) 

Tl + T 2 - 2txT 2 



(10) 



We note that the entanglement is maximum (C = 1) when 
Ti = t 2 , and minimum (C — 0) when n = and t 2 = 1 
or Ti = 1 and t 2 = 0. The scattering matrix for the 
quantum dot of Fig. [I] reads 



S 



r t' 
t r' 



(11) 



where r,r',t, and t' are 2x2 reflection and transmission 
matrices, respectively. Thus, ti and t 2 in Eq. (|10[) are 
the eigenvalues of the product ttf . In the presence of 
TRI, S is unitary and symmetric. If TRI is broken (due 
to, e.g., the application of a magnetic flux), then S is 
only unitary. 

On the other hand, the chaotic scattering in the quan- 
tum dot gives a stochastic character to the entanglement 
production. Therefore, a statistical analysis of the entan- 
glement is required. Previous works focused only on the 
mean value and variance of the concurrence!^ As we 
have mentioned, here we obtain the complete distribu- 
tion of the concurrence and entanglement of formation. 
In the framework of random matrix theory, the statistical 
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FIG. 2: Distribution of concurrence in the presence {(3 = 1) 
and absence ((3 = 2) of time reversal invariance. Insets: Nu- 
merical verification of the Pp(£). A good agreement between 
theory and numerics can be observed in both symmetry cases. 



properties of the S matrix depend on the symmetry class 
/3 i 14 i 15 In the presence of TRI, the statistics of an en- 
semble of unitary and symmetric S matrices is described 
by the so-called circular orthogonal ensemble ((3 = 1). 
When TRI is absent, the statistical properties of S are 
described by the Circular Unitary Ensemble ((3 = 2). 
Particular attention has been given in the past to the 
statistical properties of the transmission eigenvalues T n 
due to their relevance in quantum transport phenomena. 
In fact, the joint distribution P/3({T n }) is known. Here 
we are interested in the joint distribution of only two 
eigenvalues, which is given b y 16 ' 17 



P/3(n,r 2 ) = cp\ri - T 2 f(T 1 T 2 ) f3/ 



2-1 



(12) 



where cp is a normalization factor. Thus, the production 
of maximally or minimally entangled states is eventually 
determined by the statistical distributions of n and r 2 
through the relations ((HJ) and lflQ|), 

We can now calculate the distribution of the concur- 
rence Pp{C) which is defined as 



2 v V 1 (1-ti)t 2 (1-t 2 ) 



Ti + r 2 — 2tit 2 



(13) 







where (• • • ) stands for the ensemble average performed 
with the probability distribution of Eq. (frSJ ■ The double 
integral over the variables t\ and t 2 in Eq. (|13[) can be 
performed exactly. After somewhat lengthy but direct 
algebra, we obtain a surprisingly simple result for the 
case of preserved TRI {(3 = 1): 



Pi{C) 



{1 + Cf 



(14) 



whereas for broken TRI ((3 = 2), we obtain the more 
cumbersome result: 



W) = 



2C 



(1-C 2 ) 3 



3(2 + 3C 2 )arctanhVl - C 2 



(ll + 4C 2 )Vl-C 2 



(15) 



In Fig. [2l we plot the distributions P\{C) (solid line) 
and P 2 (C) (dashed line) as given by Eqs. (fT4| and lfL5f . 
respectively. The resulting curves look quite different. 
The probability of producing maximally entangled states 
(C = 1) is finite for (3 — 1, while for (3 — 2 such probabil- 
ity vanishes. Besides, separable states (C = 0) are pro- 
duced with a maximum probability for (3 = 1, whereas for 
(3 = 2 the distribution goes to zero for such disentangled 
states. 

As we noted previously, maximum entanglement (C = 
1) requires t\ = t 2 . Due to the "repulsion" factor |n — 
T2f in the joint distribution (fT2"| . one might expect that 
Pp(C = 1) = for both (3 = 1 and 2. However, this is 
not the case for (3=1 due to an exact cancellation of 
the repulsion factor when performing the integrals over 
7i and r 2 in Eq. (fl"3"f . As a verification of our analytical 
results, we compare the distribution given by Eqs. (fT4"|) 
and (|15[) with numerical simulations of a chaotic quantum 
dot as performed in Ref. H. A good agreement between 
numerical and analytical results can be seen in the left 
and right insets of Fig. [2 

From the distributions (fT4")l and (fl~5"|) . we can calculate 
all the moments of Pp(C). In particular, for the mean 
value, we obtain 

,,_/ ln4-l« 0.3863 for (3 = 1, . . 

[L) \ 47r(2l7r/64 - 1) ps 0.3875 for (3 = 2, ( ' 

while for the variance of the concurrence we have 

f 2[1 - 2(ln 2) 2 ] ps 0.0782 for (3 = 1, 
var(C) = \ 5^(88 - 2l7r)(2l7r - 40) - 22 ps 0.0565 
[ for (3 = 2 . 

(17) 

These results are in agreement with those reported in 
Refs.0 andd. Both moments are practically independent 
of whether TRI is preserved or not. From Fig.[2l however, 
it is clear that the first two moments are not sufficient 
for a complete characterization of the distributions. In 
fact, the fluctuations given by y/vax(C) are of the order 
of (C). Seminal works 7 - based on the study of only (C)p 
and var^(C) arrived at preliminary conclusions suggest- 
ing that the breaking of TRI has no significant effects on 
the production of entanglement, which is in contrast to 
other quantum properties linked to transport such as the 
conductance. By studying the complete distribution, we 
find that this is not the case and we conclude that TRI 
has remarkable effects in the production of entanglement. 

Once we have calculated the distribution for the con- 
currence, we can easily obtain the distribution for the 
entanglement of formation Qp(£). We just need to make 
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FIG. 3: Entanglement of formation distributions in the pres- 
ence (/3 = 1), and absence of time reversal invariance (/3 = 2). 



the change of variable C — > £ in Eqs. (|T4[) and (|15|) by 
using Eqs. (J8j> and ([9]). In this case, however, we can not 
give analytic expressions for Qp{£) because of the tran- 
scendental functions in Eq. §§§ . We proceed numerically, 
instead, by calculating the concurrence as a function of 
the entanglement of formation, C(£). From the Jacobian 
of the transformation C — > £, we find that Qp(£) is re- 
lated to Pp{C) through 



Qp(£) 



l 



\n2y/\-C{Sf 



-Pp(C(S)). (18) 



C(£)atanh(^l-C(£) 2 )~ 
The distribution Qp(£) is depicted in Fig. [3j Of course, 



these curves give us the same information provided by 
Pp(C) in Fig. [21 and the interpretation of the TRI effects 
on the entanglement production is similar. However, £ 
has the advantage of quantifying the cost of preparing an 
entangled state in terms of the amount Bell pairs needed. 
As regards the mean value and variance off, we obtain 



(£) = 



0.285 for 13 
0.273 for/3 



(19) 



and 



var(£) 



0.078 for 0=1, 
0.056 for j3 = 2 



(20) 



These two moments do not reveal any significant TRI 
effect, similar to what was obtained for the concurrence. 

Summarizing, we present a statistical analysis of the 
production of orbital entanglement in quantum-chaotic 
dots by calculating the full distributions of the corre- 
sponding entanglement measures. Particularly, we found 
simple analytical expressions for the concurrence distri- 
butions in the presence and absence of TRI. The com- 
plete distributions of these measures allowed us to iden- 
tify significant TRI effects on the entanglement produc- 
tion, which are not reflected at the level of the first and 
second moments. 

We thank S. Montangero for helpful discussions. This 
work was supported by the "Ramon y Cajal" program of 
the Spanish Ministry of Education and Science. 



1 A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

2 M.A. Nielsen and I.L. Chuang, Quantum Computation and 
Quantum Information (Cambridge University Press, Cam- 
bridge, 2000). 

3 C.W.J. Beenakker, in Quantum Computers, Algorithms 
and Chaos, Proceedings of the International School of 
Physics "Enrico Fermi," Varenna, 2005, (IOS Press, Ams- 
terdam, 2006). 

4 G. Burkard, J. Phys.: Condens. Matter 19, 233202 (2007). 

5 P. Recher, E.V. Sukhorukov and D. Loss, Phys. Rev. B 
63, 165314 (2001); G.B. Lesovik, T. Martin, G. Blatter, 
Eur. Phys. J. B 24, 287 (2001); W.D. Oliver, F. Yam- 
aguchi, and Y. Yamamoto, Phys. Rev. Lett. 88, 037901 
(2002); V. Bouchiat, N. Chtchelkatchev, D. Feinberg, G.B. 
Lesovik, T. Martin, and J. Torres, Nanotechnology 14, 
77 (2003); D.S. Saraga and D. Loss, Phys. Rev. Lett. 90, 
166803 (2003). 

6 S. Bose and D. Home, Phys. Rev. Lett. 88, 050401 
(2002); C.W.J. Beenakker, C. Emary, M. Kindermann, 
and J.L. van Velsen, Phys. Rev. Lett. 91, 147901 (2003); 
P. Samuelsson, E.V. Sukhorukov, and M. Biittiker, Phys. 
Rev. Lett. 92, 026805 (2004); ibid., New. J. Phys. 7, 176 
(2005); A.I. Signal and U. Ziilicke, Appl. Phys. Lett. 87, 
102102 (2005); A.V. Lebedev, G.B. Lesovik, and G. Blat- 



ter, Phys. Rev. B 71, 045306 (2005). 

7 C.W.J. Beenakker, M. Kindermann, CM. Marcus, and A. 
Yacoby, in Fundamental Problems of Mesoscopic Physics, 
edited by I. V. Lerner, B.L. Altshuler, and Y. Gefen, 
NATO Science Series II vol. 154 (Kluwer, Dordrecht, 
2004). 

8 D. Frustaglia, S. Montangero, and R. Fazio, Phys. Rev. B 
74, 165326 (2006). 

9 An orbital entangler for electrons based on an interacting 
mechanism was first introduced in P. Samuelsson, E. V. 
Sukhorukov, and M. Biittiker, Phys. Rev. Lett. 91, 157002 

(2003) . 

L. Faoro, F. Taddei, and R. Fazio, Phys. Rev. B 69, 125326 

(2004) . 

Alternatively, the degree of entanglement could be deter- 
mined by introducing some kind of state tomography as, 
e.g., in G. Burkard and D. Loss, Phys. Rev. Lett. 91, 
087903 (2003); V. Giovannetti, D. Frustaglia, F. Taddei, 
and R. Fazio, Phys. Rev. B 74, 115315 (2006); ibid., Phys. 
Rev. B 75, 241305(R) (2007). 

2 H. M. Wiseman and J. A. Vaccaro, Phys. Rev. Lett. 91, 
097902 (2003). 

3 W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

4 C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 

5 P. A. Mello and N. Kumar, Quantum Transport in Meso- 



■5 



scopic Systems. Complexity and statistical fluctuations 

(Oxford University Press, Oxford, 2004). 

H. U. Baranger and P. A. Mello, Phys. Rev. Lett. 73 142, 

(1994). 



R. A. Jalabert, J.-L Pichard, and C. W. B. Beenakker, 
Europhys. Lett. 27, 255, (1994). 



